Matching fields of a long superconducting film 
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We obtain the vortex configurations, the matching fields and the magnetization of a supercon- 
' ducting film with a finite cross section. The applied magnetic field is normal to this cross section, 

and we use London theory to calculate many of its properties, such as the local magnetic field, the 
free energy and the induction for the mixed state. Thus previous similar theoretical works, done for 
04 ' an infinitely long superconducting film, are recovered here, in the special limit of a very long cross 

section. 



PACS numbers: 74.60.Ec, 74.60. Ge 



I. INTRODUCTION 



Oh' 
^ . 

. The superconducting properties of an infinite type-II film have been studied for more than three 
! decades. As long as 1964, it was first showed by AbrikosovS that, when an infinitely long film is 
placed in an external magnetic field, parallel to its surfaces, vortices penetrate collectively into the 
sample in the form of linear chains. He also calculated the lower critical field Hci{b) as a function 
^ . of the film thickness b and showed that for small b, it becomes increasingly more difiicult _for a 
O \ linear chain of vortices to penetrate the film. After the work of Abrikosov, many theoreticaltTEl and 
experimental wor kslH have been carried out to study the mixed state of a type-Il superconducting 
^ ' film. As the external magnetic field exceeds the value of Hci{b), the linear chains of vortices start 
>- . to penetrate into the sample, one at the time. This refiects into the magnetization, which exibits 
^ [ a succession of peaks related to the penetration of each linear chain. The transition between two 
^ vortex states, differing by a linear chain, occurs at a well defined critical magnetic field, known as 
■ the matching field. 

In this paper we determine the matching fields and the magnetization of a film of finite cross 
Q\ [ section, namely, with a finite width a and also with a finite thickness b. The applied field is always 
normal to this cross section, so that the third dimension is also taken very long here. For b very 
large we recover the results obtained by the previous authors in the limit of an infinitely long film. 
Interestingly there are two kinds of matching fields in this geometry, because vortices penetrate 



. individually here. The first kind corresponds to the critical field for the penetration of an individual 
Q ' vortex. Thus there are many of this critical field, each delimiting the transition from a A^— to a 
O (A^ + 1)~ vortex state. The second kind is associated to a transition between different chain states, 
>■ . like in the case of an infinitely long film. As vortices enter the finite film, they align themselves and 

■ form a linear chain, parallel to the longer side of the film. As the number of vortices increase inside 

■ the sample, the linear chain saturates and undergoes a transition, splitting into two new chains. 
This process continues from two to three chains and so on. The the critical fields associated to these 
transitions correspond to that required to enter a new linear chain in the case of the infinitely long 
film. 

In summary we show here several new properties of the finite film, such as in the formation of 
the equilibrium vortex state. The source of these new effects lays in the strong and very important 
border effects that do not exist in case of the infinitely long film. For instance, a vortex state made 
of parallel chains is substantially affected near the edges of the longer side of the film. In order to 
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allow for any possible configuration, we let the positions of the vortices entirely free and minimize 
the Gibbs free energy with respect to them using a numerical procedure. For this end, we use the 
Monte Carlo Simulated Annealing minimization method. This method allows us to search for the 
global minimum of the Gibbs free energy (or at least very close to the global minimum) of the vortex 
system, which corresponds to the most stable configuration. We also solve here for the matching 
fields, which requires previous knowledge of the Gibbs free energy minimimum for a fixed number 
of vortices. Notice that this program is much simpler for an infinite film, because there one can 
suppose some type of symmetry for the vortex state, and minimization can be reduced to a small 
set of parameters, such as the lattice spacing and the number of chains^. 

Our analysis is restricted to the low magnetic field regime and to the very strong type-II super- 
conductors, where the vortex-cores do not overlap and the distance over which the magnetic field 
penetrates into the sample is much larger than the size of the cores. Within this regime, the mixed 
state may be well described by the London theory accompanied by some regularization procedure 
to cure some divergencies which appear at short length scales. To the best knowledge of the present 
authors this study has not been done before. 



II. THE MIXED STATE 



The description of the mixed state of a type-II superconductor involves two fundamental lengths. 
One is the distance over which the magnetic field penetrates into the sample, the penetration length 
A. The other one is the coherence length, the size of the vortex-core. Our treatment will be limited 
to the case of a strong type-II superconductor where the Ginzburg-Landau parameter k = A/^ 1. 
Vortices interact strongly since the distance between them can be much smaller than A. Their 
repulsive interaction forbids overlapping, namely the onset of vortices with multiple vorticity. In 
this context, the distribution of the local magnetic field may be well described by the London 
equation. In rectangular coordenates (x, y) the London equation for the local magnetic field h = hz 
is given by 

- AV^/i + = $0 5: 5(r - r,) , (1) 

i 

where $0 is the quantum flux and rj = (xj, r/i) is the position of the i-th vortex. We have assumed 
that the vortices are straight lines. Therefore, a 3D problem is reduced to a 2D one. 

The geometry of the problem is ilustrated in Fig. |I|. Let a and b be the sides of the superconducting 
film. We will solve this equation subject to the following boundary conditions 

= h{x,b) = H, 

= , (2) 

y=0,b / y=±a/2 

where H is the magnitude of the external magnetic field at the interface vacuum-superconductor and 
is pointing along the z-direction. The second condition imposes that the supercurrents are confined 
into the sample. 

Other authors have used the image method to solve the London equation (See for instance Ref. 
^ and p. We prefer to solve this equation directly, by using the Green's function method. The 
equation for this function associated with the London equation is given by 

- X^V^G + G = 5{x - x')5{y - y') , (3) 
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where G must satisfy the following boundary conditions on the film borders 



G{±a/2, y) = G{x, 0) = G{x, b) = 



(4) 



We can find an expression for the local magnetic field in terms of the Green's function by multipling 
(H) by G{x,y,x',y') and @ by h{x,y), subtract the results and then integrate what is left. One 
obtains 



h{x\y') = 00 XI G{xi,yi,x\y') 



-H 



/-a/2 rb 

I- dx dyG{x,y,x\y') 

J-a/2 JO 



(5) 



where we have made use of the boundary conditions @ and (H). 

Then, the solution for the local magnetic field is transfered to the determination of the Green's 
function. The approach used to find the solution for this function is sketched in full details in Ref. 



11, except by the fact that there the boundary conditions are taken at infinite. Here we only show 



the main steps of how to calculate this function. First of all we expand G in a Fourier series 



G[x, y^x,y) = - 2^ sm I — — 

m=l \ / 



sm 



rmry 



(6) 



which satisfies the boundary conditions of (^) at y = 0, b. 
Inserting this equation into (^) we find for gm 



ym I ^2 



dx"^ 



5{x — x') 



(7) 



where we have used the fact that the sequence {yf sin (^^^^ 
orthonormal functions, that is. 



m 



1, 2, 3 . . .| is a complete set of 



2 ^ . (niTxy'" 

5 S I—, 



sm 



rm^y 



Ky - y') ■ 



Here 



.1/2 



(9) 



The function gm{x,x') must satisfy the same boundary conditions as G{x,y,x',y'), that is, 
(7m(±'^/2, x') = 0. In addition, the derivative of gm.{x,x') is discontinous at x = x' which results 
from the delta function on the right-hand side of (|^). 

Under these conditions, the solution for gm{x,x') is given by 



9m{-^i X ) 



2Xam sinh(ama/ A) 



{cosh[a;m(|x — x'\ — a)/A] — cosh[am(x + x')/A]} 



(10) 



This completes the solution for the Green's function. To proceed we still have to find the integral 
of the Green's function which appears in (||). This involves a tedious, but straightforward algebra. 
One has 
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j-a/2 j-b 4 °o 

/ dx dyG{x,y,x\y') = - ^ 

J-a/2 Jo ±^ 



:-iy 



sm 



X 



cosh(a2m+ia:VA) 
cosh(a2m+io./2A) 



Notice that this vanishes at x' = ±a/2 or y' = 0,6, so that the second term of (H) is just H as 
required by the boundary conditions of @. 



Next we can simphfy (|TT]) by noticing that 
'-D™ b 



A 

4 y 

b aL+i (2m + 1)77 



sm 



rriTcy 



dygiy,y') 

cosh[(y'-6/2)/A] 
cosh(6/2A) 



(12) 



where 



, 2 ~ 1 . (mwy'' 



sm 



/miry 



= 2Asinh(fe/A) {^Qsh[(||/ - ^)/A] - cosh[(|/ + yO/A]} • 
By combining (^), ([TTD and (^2|) we finally obtain for the local magnetic field 

/ cosh[(?/ - b/2)/X] 



(13) 



h{x, y) = ^oJ2 Gix, y, Xj, yi) + H 

i 

4 ^ f-l)'" b 



V 



^ & £) "L+i (2m + l)7r 



sm 



cosh(6/2A) 
(2m + l)7r?/ 



cosh(a2r»+ia^/A) \ 
cosh(a2m+ia/2A) J 



(14) 



where we have made use of the following symmetry property of the Green's function, G{x, y, x', y') = 
G{x', y', X, y). We can see that the local magnetic field is composed essentially by three contributions. 
The last two terms, proportional to H in (|1^), represent the penetration of the external magnetic 
field near the surface when the vortices are still absent. Once the vortices start penetrating the 
sample they will be trapped by the shielding currents associated with the last two terms of (0). 
The Green's function in (|1^) contains two terms (see (^) and (jl^)). The first one, which depends of 



represents the local magnetic field of the vortices located inside the sample, and the second 
one, which depends on x + is the local magnetic field of the image vortices located outside the 
sample. 

The free energy per unit volume of the superconducting film is given by 



1 



SttA 



+ 



a/2 
-a/2 




x=a/2 



'dh' 

dx 



-a/2^ 



dx -— 



(15) 



where A = ab. 
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Now notice that by integrating the kinetic term of the London equation of (|I]) we obtain precisely 
the second term of (p^Sf). Then, this equation can be simphfied to 



(16) 



where N is the number of vortices, B is the spatial average of the local magnetic field. The induction 
is obtained by integrating (|I^) over the cross section of the superconducting film. One has 

cosh(?/i - 6/2)/A) 



AB = iV$r 



4$o 



i 

oo 

EE 



cosh(6/2A) 
-D™ 1 



+HA 



vr a^+i (2m + 1) 

J tanh(6/2A) 8 



sm 



(2m + l)7ryi 



E 

m=0 



(2m + l)a2m+i 



COsh{a2m+lXi/X) 

cosh(a2m+ia/2A) 

tanh(a2m+ia/2A) 1 
(a2m+ia/2A) J 



(17) 



Upon substituting (0) into (|I6|) we are left with 



SttA 



J2 G{xi,yi,Xj,yj) 



+ 



J tanh(6/2A) 8 
8^ 1 (6/2A) ^ 



E 

m=0 



H 2 



(2m + 1)^2 



m+l 



tanh(a2m+io/2A) 
(a2m+ia/2A) 



In order to study the most stable configuration of the vortex lattice, it is more convenient to take 
the Gibbs free energy. In units of volume, its expression is given by ^ = JF — BH/Att. From (|lj 
and (|T8|), it follows that 



87iA 

>0 



■kA 



EE 



AtiA Y cosh(6/2A) 
1 . \{2m+l)'Kyi 



«2m+i (2m + 1) 



sm 



[ tanh(6/2A) _ _8_ ^ 
'8^ I (6/2A) ^ ^ 



m=0 



(2m + l)a2m+i 



cosh(a2m+ia^i/A) 
cosh(a2m+ia/2A) 

tanh(a2rrt+ia./2A) 
(a2m+ia/2A) 



-N 



AttA 



(19) 



The last two terms are the energy of the Meissner state as if no vortices were present. The second 
and third contributions are the potential barrier which pin the vortices inside the sample. The 
Green's function in the first term gives rise to two contributions. One is the repulsive interaction 
between the vortices and the other one is the atractive interaction between the vortices and the 
images which are virtually placed outside the sample. 

We shall work in the limit of thin film in which (vrA/fe)^ ^ 1. Within this limit, the sum in m can 
be evaluated exactly. We proceed as in Ref. H. One has. 



Cr(^Xi, yi, Xj, yj) 



47rA2 



In 



cosh[7r 


-x,\/h] 


- cos[7r(?/i + yj)/h] 


cosh[7r Xj 


-xm 


- cos[7r(|/i - yj)/h] 
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In 



In 



cosh[7r(a — Xi — Xj)/h] 



cos 



+ %■)/&] 



cosh[7r(a 



cos 



cosh[7r(a + Xi + Xj)/h] — cos[7r(yj + yj)/h] 
cosh[7r(a + Xi + Xj) /h] — cos[7r(yj — yj)/h] 



(20) 



The first term in ( ]20| ) is the same as the one found in Ref. |^. The other two terms are due to the 
presence of the lateral border of the film. London theory is not regular for vortex self-interaction. 
In fact, one can notice that the contribution i = j to the Gibbs free energy (|TUp gives rise to a 
logarithmic divergence. We remediate this divergence by using a sharp cutoff in which | X Xj I IS 
replaced by ^ for i = j. In the limit (6/^0 ^ 1 we obtain 



G{xi,yi 



47rA2 
-In 



In 



cosh[7r(a 



(vre/6)2 
2xi)/b] - 



cos{2nyi/b) 



In 



cosh[7r(a — 2xi)/b] — cos(7r^/6) 
cosh[7r(a + 2xi)/b] — cos(27r?/j/6) 
cosh[7r(a + 2xi)/b] — cos{n^/b) 



(21) 



The sum over m in the third term on the right-hand side of (|T8D cannot be evaluated exatly in the 
limit {nX/b)^ ^1. So it will be kept as it stands. 



III. RESULTS AND DISCUSSION 



The lower critical field, the magnitude of the external field sufficient to create a vortex inside the 
sample, may be obtained by neglecting the quadratic term in (^) and equating the remaining terms 
to zero. We find 



Hci{a,b) 



$0 

2A2 



^lnf^)-ln 



cosh(7ra/fe)+l 
cosh(7ra/fe) — 1 



1 



1 



cosh{6/2A) 



-4E 



(-1)'" 1 

m=0 (2m+l)7ra^^_^^ cosh(a2m+ia/2A) 



(22) 



where we have assumed that a single vortex is at the center of the film. It can be easily seen that 
for a ^ cxD we obtain the same result as in Ref. ^ A quick inspection of ( |2^ ) shows us that as 
a and/or b decreases, Hci{a,b) increases. In other words, size effects provoke a delay in the first 
penetration of flux, assuming the external field increasing with time. These finite size effects have 
been observed in supercondutors of different geometryGj. The determination of the other critical 
fields will be discussed below. 

Now we turn our discussion to the minimization of the Gibbs free energy with respect to the vortex 
positions. Our investigation will be guided by the Monte Carlo Simulated Annealing minimization 
method. We have done this from N = 1 up to = 30. We start from an arbitrary initial 
configuration of the vortex positions at a certain "temperature". By adiabatically lowering this 
ficticious temperature. Simulated Annealing will concentrate on the region where it is most likely to 
succeed. 

The parameters we use here are k = 100, a = 7A, and b = 0.75A. We find that, as the penetration 
of flux is initiated, the vortices will be symmetrically located with respect to the center of the film in 
form of a linear chain aligned parallel to the x-axis (the longer side of the film). As H is increased, 
at a certain critical value of the external magnetic field, this linear chain breaks into a double chain. 
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At the interior of the film, the vortices arrange in form of a "zig-zag" . However, at the edges of the 
film both linear chains are distorted and joined together (see Fig. 0). For the parameters quoted 
above, we have found that the change in the symmetry of the lattice occurs at = 24. This critical 
value of is expected to grow in case the value of the width a of the film becomes larger. 

In Fig. ^, we also have shown some other configurations for higher values of A^. Notice that the 
vortex state develops some unusual patterns in which the linear chains entangle in the interior of the 
film. On each step of our numerical calculation, to search the minimum G we lowered the temperature 
as slowly as possible. However, we cannot assure that the minimum G found corresponds to the 
global optimum. So, those entangled chain configurations may be metastable solutions. 

The penetration of each additional vortex corresponds to a well defined value of external magnetic 
field which we denote by Hg]^; H^i = Hci{a, h). We find these matching fields by assuming that at the 
transition from A^ to A^+1 vortices, the Gibbs free energy is continuous, that is, Qjsi = Gn+i- Because 
the vortex positions depend of the external magnetic field, this is a trancendental equation. Full 



details of how to solve this transcendental equation can be found in Ref. |12[ We have determined the 
sequence of matching fields for the parameters quoted above. We then used these values to calculate 
the magnetization, which is defined by M = {B — if)/47r. 

The magnetization as a function of H is depicted in Fig. ^. At the transitions, the creation of 
another vortex will force a rearrangement in the chain. This in turn will provoke a discontinuity of 
the magnetization at if = Hsn- We can see that the magnetization changes in small steps. So, the 
magnetization can be used to measure the number of individual vortex penetration in the sample. 
Each peak indicates that a new vortex entered the film. 

Notice that —AnM has a maximum for H > Hci{a,b), signaling a transition from a single to 
a double chain state, in constrast to superconducting systems with no edges, that just show a 
monotonic behavior. In each phase between if^Af and ifg^Ar+i, for each value of H we should calculate 
the corresponding equilibrium configuration of the vortices and use it to evaluate the induction B. 
For the parameters used in Fig. ^ the differences AH = Hg^N+i — Hsn are not to large. Then, the 
vortex positons should not vary significantly with H within each phase. So, to calculate B, we held 
the vortex positions fixed at if = Hsn- As a consequence, the magnetization varies linearly with H 
within each phase. 

We can also observe that the intensities of the peak of the magnetization at the transition from a 
single to a double chain lattice, and aX H = Hd {a, b) have approximately the same intensity. In Fig. 
§ we repeated the same calculation with a = 7X and b = 0.5A. For A^ up to 30 vortices we did not 
find any break in the linear chain (see Fig. However, the magnetization curve shows a tendency 
to reach the second peak of slower magnitude than at if = Hd- On the other hand, for an infinite 
film, as shown in Ref. ^, the intensity of the magnetization at the phase transitions tends to grow as 
H increases. We then suspect that this behavior of the magnetization may be a result of size effects. 

Finally, we would like to point out that the successive discontinuities of the magnetization discussed 
here should be experimentally observed. To see this, we note that the order of magnitude of AH = 
ifs^AT+i — HsN in Fig. |^ can be as high as $o/2A^. For zero temperature, A is typically of order 10^ 
Awhich gives Aii ~ 10^^ T. On the other hand, the jump in the magnetization can achieve the 
magnitude of order 10^'' T. This is well inside the resolution of a magnetometer. 
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FIG. 1. Cross section of the superconducting film. The external magnetic field is along the ^;-axis and its value is H for 
\x\ > a/2, and y < and y > b. 

FIG. 2. The double chains for A'^ = 24, 25, 26, 27, in the clockwise sense. The parameters used are quoted in the text. Notice 
that the inner vortices make a "zig-zag" as in a triangular lattice. However, at the edges of the film, both chains end up at the 
same vortex. 



FIG. 3. The magnetization as a function of the external magnetic field H. The parameters used are quoted in the text. 



FIG. 4. The magnetization as a function of the external magnetic field H. The parameters used are quoted in the text. 
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